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Abstract 



p ■ We construct two new exactly solvable potentials giving rise to bound-state so- 

lutions to the Schrodinger equation, which can be written in terms of the recently 
Ch ' introduced Laguerre- or Jacobi-type Xi exceptional orthogonal polynomials. These 

^ ■ potentials, extending either the radial oscillator or the Scarf I potential by the ad- 

O^i dition of some rational terms, turn out to be translationally shape invariant as their 

standard counterparts and isospectral to them. 
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Classical orthogonal polynomials are known to play a fundamental role in the construc- 
tion of bound-state solutions to exactly solvable potentials in quantum mechanics. For such 
a purpose, the factorization method Pfl E] and its realization in supersymmetric quantum 
mechanics (SUSYQM) [3], especially for shape-invariant potentials [1], as well as the equiv- 
alent Darboux transformation [5j , prove very useful. The same is true for a more traditional 
approach, the point canonical transformation (PCT) method [6j, consisting in directly map- 
ping Schrodinger equations into the second-order differential equations satisfied by those 
polynomials. 

On the other hand, bound-state solutions to exactly solvable potentials are by no way 
restricted to classical orthogonal polynomials. For instance, SUSYQM [3 [H El |10] and 
the Darboux transformation [Til [l2l [IS] are very efficient at producing new sophisticated 
exactly solvable potentials by adding or deleting some states, or else by leaving the spectrum 
unchanged. The PCT method is also very powerful for generating new shape-invariant or 
non-shape-invariant potentials not only in a standard context [15], but also in more general 
ones, such as those of quasi-exact p£j or conditionally-exact p/Tj solvability and that of 
position-dependent masses [18\ . 

Very recently, two new families of exceptional orthogonal nth-degree polynomials, 
P^'^\x) and L!n\x), n = 1, 2, 3, . . . , have been introduced [191 [20]. Such sequences, 
referred to as Jacobi- or Laguerre-type Xi polynomials, respectively, arise as solutions of 
second-order eigenvalue equations with rational coefficients. They are characterized by the 
remarkable property that although they do not start with a constant but with a linear poly- 
nomial, they form complete sets with respect to some positive-definite measure in contrast 
with what would happen if one deleted their first member from the families of classical 
orthogonal polynomials. 

In this communication, we plan to show that there exist some exactly solvable potentials 
whose bound-state wavefunctions can be written in terms of these new exceptional orthog- 
onal polynomials. For such a purpose, we shall make use of the standard PCT method. In 
a second step, we shall employ SUSYQM techniques to prove that our new exactly solvable 
potentials are transitionally shape invariant. 
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In the PCT method, one looks for solutions of the Schrodinger equation, 

Hijix) = (^-^ + Vix)^ ^Ijix) = E^Pix), (1) 

of the form 

^(x) = fix)Figix)), (2) 

where f{x), g{x) are two so far undetermined functions and F{g) satisfies a second-order 
differential equation 

F + Q{g)F + R{g)F = 0. (3) 

Here a dot denotes derivative with respect to g. 

On inserting equation ([2]) in equation ([1]) and comparing the result with equation Q, 
one arrives at two expressions for Q{g{x)) and R{g{x)) in terms of E — V{x) and of f{x), 
g{x) and their derivatives. The former allows one to calculate /(x), which is given by 

1 / 1 /"^("^^ \ 
f{x) oc ^ exp - / Q{u) du , (4) 



V \2 

while the latter leads to the equation 

2 

2g' 4 \g') ' ^ 2^ 4 

In (jlj) and ([5]), a prime denotes derivative with respect to x. For equation (JSj) to be satisfied, 
one needs to find some function g{x) ensuring the presence of a constant term on its right- 
hand side to compensate E on its left-hand one, while giving rise to a potential V{x) with 
well-behaved wavefunctions. 

Let us start by considering for the second-order differential equation ([3]) that satisfied 
by Laguerre-type Xi polynomials L^n\x), n = 1, 2, 3, . . . , a > 0. In such a case, the 
functions Q{g) and R{g) can be expressed as [191 EO] 

{g - a){g + a + 1) a + 1 2 

Q{9) = 7 — : — ^ = -1 + 



g{g + a) g g + a 

1 fg-a \ n-2 2 

Rid) = - + n - 1 = + , 

g \g + a J g g + a 
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so that we obtain 

1. _1 2__1 2an + a'^-a + 2 1 + 2 

2^ 4^ ~ i+ ^^^^ + V (^T^' ^ ^ 

A constant term can be generated on the right-hand side of equation ([5]) by assuming 
g''^/g = C, which can be achieved by taking g{x) = jCx^. Equations ([5]) and ([6]) then yield 

E = ^C{2n + a-l), 



On setting 



we arrive at 



and 



16 x2 Cx^ + Aa {Cx^ + Aaf 



C = 2iu, a = I + ^, n = u + 1, 



= cu(2z/ + / + |), z/ = 0,1,2,..., (7) 



V{x) =Vi{x) + V2ix), 

Vi{x) = -Uj\^ + 



V2{X) 



4 x^ 

Auj 8uj{2l 



ujx'^ + 21 + 1 + 21 + If 



For 



< X < oo, w > 0, / = 0,1,2,..., 

V{x) is a well-behaved potential, which may be interpreted as an /-dependent (effective) 
potential, extending the standard radial oscillator potential Vi(x) by the addition of some 
rational terms. Such terms do not change the behaviour of the conventional potential for 
large values of x, while for small values they have only a drastic effect when the angular 
momentum / vanishes, in which case V{Q) = —Au < instead of V{0) = 0. 

From equation ([7]), it follows that the extended potential has the same spectrum as the 
standard one. The corresponding wavefunctions can be found from equations ([2]) and (jl]). 
On solving the latter for the choices made here for Q{g) and g{x), we get 
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where the normahzation constant is obtained from equations (31), (33) and (34) of [19] as 



, \ 1/2 



In particular, the ground-state wavefunction can be written as 

1 2 
oc i)w{x)[l + 0(x)], oc x'+^e"^"^^', 0(x) = ^ 2/ + 1 ' 

and differs from that of the standard radial oscillator, ^/'io(x), by the extra factor 1 + 0(a;). 
It is obvious that it is a zero-node function on the half-line, as it shoud be. Furthermore, 
it can be easily checked by direct calculation that it satisfies equation ([T]) for the potential 
dHD and Eo = cu (/ + §). 

More generally, as shown in [20], the polynomial L^j^-^^ (J^ux^) (and hence the wavefunc- 
tion il)y{x)) has V zeroes on the half-line. From general properties of the one-dimensional 
Schrodinger equation, it therefore results that we have found all the eigenvalues of potential 



Let us next consider the case where the second-order differential equation ([3]) coincides 
with that satisfied by Jacobi-type Xi polynomials P!C'^\x), n = 1, 2, 3, . . . , a, /3 > —1, 
"7^/3 [laEni, i.e., 

{(3 + a + 2)g-{[3-a) 2{p - a) 



Qig) 

m 



l-g' iP-a)g-{P + ay 

{l3-a)g-{n-l){n + (3 + a) - a)^ 



{(3-a)g-{(3 + ay 
This choice leads to 



2" 4" l-g2 (I-9T (13 - a)g - (P + a) [(/3 - 0)9 - (;3 + a)P ' 

where 

G=^{f3-a){f3 + a), J = -^{[3^ + a' - 2) , 
K= , L = -2(/3-a). 



We can obtain a constant term on the right-hand side of ([5]) by assuming g''^ / [l—g"^) = C. 
For (7 = > 0, we can take g{x) = sin(ax). On rescahng the variable x, the parameter a 
can be set equal to one. Then with the changes of parameters and of quantum number 

a = A-B-l, [3 = A + B-\ or A=|(/5 + a + l), B = \{(3-a), 
n = u + 1, 

we arrive at the following results 

E, = {u + Af, u = 0,1,2,..., (10) 



and 



Vix) = Viix)+V2ix) 



Vi{x) = [A{A - 1) + B"^] sec^x - B{2A - 1) secxtanx, (^1) 
y(._ 2(2A-1) 2[i2A-ir-AB'] 
' 2A-l-2Bsmx {2A-l-2Bsmxy' 

The function Vi (x) defines a Scarf I potential, for which it is customary to assume 

TT TT 

— <x < -, < B < A-1. 
2 2' 

For such values of the variable and of the parameters, the full potential V{x) has the same 
behaviour as Vi{x) for x ±7r/2, only the position and the value of the minimum being 
modified. 

It results from ffTOl) that on using as usual Dirichlet boundary conditions at the end 
points of the interval, the extended Scarf I potential (fTTj) has the same spectrum as the 
conventional one. Its wavefunctions can be written as 



(1 - sinx)5(^~'^)(l + sinx)5(^+^) ^(a-b-^.a+b- 
2A-1 - 2Bsmx 



1 ' 



'4)y{x) = My — ^ ^ „ . P^Vi ^' (sinx) 



where the normalization constant 

.r_Bf v\{2v + 2A)V{v + 2A) ^^'^ 

"~ 2^ \ {u + A - B + \) {u + A + B + \)V {p + A - B - \)V {p + A + B - \) 

is a consequence of equations (23), (25) and (26) of [19]. 



The ground-state wavefunction assumes the simple form 

ifjoix) (X ipio{x)[l + (f){x)], ipio{x) oc (1 - sinx)5(^-^)(l + sinx)5(^+^), 

^( ^ - 2 (12) 

'^^^^ 2A-l-2Bsmx' 

where the presence of is due to the rational terms in V2{x). As in the previous case, 

the function ipoix) has no node on the interval of variation of x and can be easily checked 

to satisfy equation ([1]). 

More generally, the polynomial -^v+i ^' ^ (sinx) (and hence the wavefunction 

i^v{x)) has u zeroes on (— f , f ) [20J, so that no other eigenvalues than flTUl) may exist for 

potential (ITTi) . 

Let us finally combine our results with SUSYQM methods [H [H [9l [10]. In the minimal 
version of SUSY, the supercharges Q and Q'^ are generally assumed to be represented by 
Q = Aa^, = where a± are combinations a± = (Ji±i(j2 of the Pauli matrices and A, 



A^ are taken to be first-derivative differential operators, A = ■^_ + W{x), = + 



x] 



with W{x) known as the superpotential. The supersymmetric Hamiltonian Hs = {Q, Q^} is 
diagonal, i.e., Hg = diag(iJ^"'"), H^~'>), and its components H^"^^ can be written in factorized 
form in terms of A and A^, 

= A^A = — ^ + V^+\x) - E, = AA^ = + V^-\x) - E, 

dx^ dx'' 

at some arbitrary factorization energy E. The partner potentials V^'^\x) are related to 
W{x) through = W\x) T W'{x) + E. 

The spectrum of is doubly degenerate except possibly for the ground state. In the ex- 
act SUSY case to be considered here, the ground state at vanishing energy is nondegenerate. 
In the present notational set-up, it belongs to H'^'^\ 

H^+^iP^'^\x) = 0, ^t\x) oc exp ( - [ W{t)dt] . 



Let us identify V^~^\x) with either the extended radial oscillator potential ([8]) or the 
extended Scarf I potential ffTTj) and take for the factorization energy Eq = uj{l + ^) or 
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Eq = A^, respectively. Then tljl^~^\x) is given by equation or equation f[T^ . The 
corresponding superpotential W{x) = —ipQ^x) / ipoi^x) can be separated into two parts, 



W{x) = Wi{x) + W2ix), Wi{x) = W2ix 



^lo' ' ' 1 + 

where Wi{x) is the superpotential for the conventional potential, i.e., 

Wi(x) = -cox — ^ ~^ or lyi(x) = Atanx — -B secx, 
2 X 

and the additional term W2{x) can be written as 

Wo(x) = 2ujx' ^ ^ 



cux"^ + 21 + 1 ujx^ + 2l + 3 

or 

/ 1 1 

W2(x) = -25 COS X — 

\2A-l-2Bsmx 2A + 1 - 25 sinx . 
respectively. 

From this, it follows that the partner potential V^~\x) to V^'^\x) (with one less eigen- 
value) is given by 

F(^)(a;) = V'^+\x)+2W'{x) = Vt\x) + Vt\x), Vl'\x) = v}^\x) + 2W^{x), i = 1,2. 

As it is well known, v} \x) is a standard radial oscillator (resp. Scarf I) potential with 
I replaced by / + 1 (resp. A replaced by A + 1). It is straightforward to convince oneself 
that a similar property relates V2~\x) with V2^\x). We therefore conclude that the two 
extended potentials ([8]) and (ITT!) are translationally shape invariant as their conventional 
counterparts. 

Summarizing, we have constructed some exactly solvable potentials for which the re- 
cently introduced Laguerre- or Jacobi-type Xi exceptional orthogonal polynomials play a 
fundamental role. Furthermore, we have demonstrated that these new potentials are shape 
invariant. It is rather obvious that the method described here could be used for other 
choices of function g{x) in order to generate other types of potentials connected with such 
polynomials. Another interesting open question for future work would be the origin of the 
(strict) isospectrality observed between the extended and conventional potentials. 
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